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Abstract 

The recent progress in the study of Gahleons, i.e. equations of sec- 
ond order with an action invariant under a Galilean transformation is 
related to work on 'Universal Field Equations' [T3] which are second or- 
der equations arising by an iterative procedure from arbitrary Lagrangians 
of weight one in their first derivatives. It is pointed out that the Galileon 
is simply a Kaluza-Klein reduction of a Universal Field Equation. An 
implicit solution to the equation of motion is presented, and a class of 
explicit solutions pointed out. The multi-field extensions of both types 
of equations are derived from a first order formalism, which is simply the 
substantive derivative of fluid dynamics. 



1 Introduction 

Recently there has been a spate of articles on the subject of Galileons, i.e. 
equations, second order in derivatives of a field tt with an action invariant under 
the infinitesimal transformation 

tt{x) H> tt{x) + c + v^j^x^ + 7r(a;)u''9^7r(x) (1) 

with x'^ space-time co-ordinates, a constant vector and c a constant. [I] 
[5] [3]. In the literature two versions of the Lagrangian in 3+1 dimensions are 
proposed; one starting with the Lagrangian 



C2 = -v/TTW, (2) 

which is invariant up to a divergence under the transformation ([T]) and con- 
tinuing with a heierarchy of increasingly complicated Lagrangians [2j, ending 
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with 



-2[[n3]-3[n][n]2)--7^(4[n][^4] 



(3) 

See also [4]. The other starts with the Klein Gordon Lagrangian |5]. Here the 
notation H^^ = df^d^ir; [H]"] = trace(n") and [tt"] = dn.W^-^ .dir; i.e [tt^] = 
d^ird^d^'Kdu'K. The factor 7 is the inverse of 1/1 + (97r)2. The resulting equation 
of motion is £5 = where 



^5 = 7' ([n]^ 



6[n]2[n2] 



8[n] [n^ 



6[n^ 



(4) 



We shall discuss this procedure in the light of previous work on iterated La- 
grangians, and relate it to these recent studies. This equation of motion can be 
very much more simply expressed thanks to a result on determinants, which is 
published in Lowell Brown's book, [6| but was independently found by myself; 
here illustrated for 4x4 matrices. Consider the determinant 
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(5) 



where T„ = trace (T"), and T is a matrix of dimension 4x4. The expansion of 
this determinant is just 4! times the characteristic polynomial of the matrix T, 
and the coefficient independent of x is just the determinant of T. 



det(T) = -6 * T4 + 8 * Ts * Ti 



3 * - 6 * * 



But, if Tui, = 



n„ 



this is just the operative factor in £5, so the 



equation of motion simply obtained by setting the Hessian det jH^^ | = 0. This is 
also called the homogeneous Monge- Ampere equation and a general solution has 
been given by A.N. Leznov and myself f^J. I recently recognised an equivalent, 
and perhaps more elegant solution in a book of T.H. Chaundy, quite some time 
ago [10]. 



1.1 Chaundy's solution 

The principle is illustrated by the 3x3 case. To solve 



det 



f^xy 



i^xy 
t>yz 



0, 



(6) 



choose four arbitrary functions /(u, w), g{u,v), h{u,v), u) constrained by 
the three relations 

xf{u, v) + yg{u, v) + zh{u, v) + k{u, v) = (l){x, y, z) (7) 
xj{u, v)u + vgu{u, v) + zhu{u, v) + ku{u, v) = Q (8) 
xfv{u, v) + ygv{u, v) + zhy{u, v) + ky{u, v) = 0. (9) 
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Here subscripts denote partial differentiation with respect to u, v. Then the 
imphcit solution of these equations for (j){x,y,z) is a solution to ([6]). This is 
easy to see; the three equations imply 

d(f>{x,y,z) d(l){x,y,z) d<j){x,y,z) 

— — ' ^^('''") = — — ' = — q; — ■ (10) 

Since the right hand sides are three functions of two variables, there must be 
a relationship amongst them; i.e. there is a functional relationship among all 
the first derivatives of say F{(j)x, (py, (I'z)- As we know this means that (j) 
satisfies the Monge- Ampere equation. This is easily seen from the assumption 
of a relationship of the form 

=0 

Differentiation with respect to the variables x,y,z respectively gives the three 
relations 

dF ^ dF ^ dF ^ „ 



^•rx '-'Vy '^^z 

dF ^ dF ^ OF ^ 

-^(Pxy + -^4>VV + TT—'Pzy 

a(j)x a(j)y 
dF , OF , dF 



(12) 



Wx''''' ^ Wy""'' ^ Wx^'' - ^^^^ 

The eliminant of these equations is simply the Monge- Ampere equation ([6|) . We 
can also verify directly this by constructing the determinant of the matrix of 
second derivatives 



det 



fu^x fv'^x fu^y ~t~ fv'^y fu^z ~^ fv^z 
guUx + QvVx guUy + 9vVy OuUz + OvVz 



(14) 



Evaluation of this determinant gives zero. It is also easy to see, by a similar 
technique that any homogeneous function of weight one, i.e which satisfies 

90 d(j) d(j) 
ox oy oz 

will also satisfy the equation (jB]) by a similar technique of differentiation by x, y, z 
in turn. Interestingly, it is easy to see in this case that the Galilean transfor- 
mation (IlJ, not now necessarily infinitesimal also satisfies the equation (fT4)) as, 
apart from the irrelevant constant, the additional terms are also of weight one, 
if TT is. One might worry as to whether the equality of the two different ways of 
evaluating mixed second derivatives introduces further constraints; for example 

5^0 _ df{u,v)du ^ df{u,v)dv ^^^^ 



dx, dy du dy dv dy 

|9^0 ^ dg{u,v) du ^ dg{u,v) dv 
dy,dx du dx dv dx 
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It turns out that the equahty of the two expressions is automatically satisfied, 

and goes back to the equality of the mixed derivatives of f,g, etc. Putting 

it more simply fuUy + f^Vy — guUx + gvVx- This is then a solution to the 
3-dimensional Galileon equation 



det 



0, 



(17) 



since this is just 6 times the determinant (jl4l) in the case where fj,, v run 
from 1. . .3. If the range is 1. . . n, the determinant decomposes into the sum of 

Hessian determinants, and the same form of solution will hold by the 

elimination of u, v from the equations 



/i— 1.. .n+1 



E 

/J.— l...n+l 



(18) 
(19) 
(20) 



The construction in general can be easily inferred, and the analogous treatment 
of the solution of a related equation, described in the following section has been 
already described [12^ . 



2 Iteration of Lagrangians 

This version of a Galileon equation is actually simply a Kaluza Klein reduction 
of what J.Govaerts, A. Morozov and I called a 'Universal Field Equation' [T3] . 
We called it this because, starting with any Lagrangian dependent only upon tt 
and its first derivatives, and homogeneous of weight one in these derivatives, i.e. 

djC 



E 



subject to an iterative procedure, gives a UFE which is independent of the form 
of the initial Lagrangian. This obtained by iterating the Euler operator £ acting 
on the Lagrangian 

d d d 
^ +^^^ ^^^JTr~ 

OTT OTTi OTTi 



£ = -—+^^—~^,^,— ... (21) 



(In principle the expansion continues indefinitely but it is sufficient for our 
purposes to terminate at the stage of second derivatives (jjij, since it turns out 
that the iterations do not introduce any derivatives higher than the second). We 

may start with a Lagrangian say, £ = / {d^n)^. Then the 5 fold iteration 

Y Ai=1...5 

£C£C£C, ■■■,££ (22) 
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where each Euler operator acts on everything to the right yields the Universal 
Field Equation 
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(23) 



This equation possesses some interesting invariance properties; First of all it is 
signature blind; being the same in a Euclidean space-as a Lorentzian space-time. 
Also, any function of a solution is also a solution! This may be verified directly, 
but it arises from the solution procedure explained in [12]. Now in a Kaluza 

di: 

Klein reduction we may choose a gauge in which — — = 1 . In this case the last 

row and column of (|23p have only two nonzero elements in the (1,6) and (6,1) 
positions, and these are unity. The resulting determinant simply reproduces 
([5]), the Hessian. 



3 Connection with earlier results 

These results are easily understood as examples of the main results of [T3] 
quoted here; 'In two recent papers |131I15| . the first with A. Morozov, hierarchies of 

Lagrangian field theories with the following properties were introduced. 

1. i) In any of these hierarchies, the Lagrangian at any given level - except of 
course at the first level - is essentially proportional to the equations of motion 
of the Lagrangian at the previous level (hence the name Euler hierarchies). 

2. ii) The proportionality factor mentioned in i) is essentially the very first La- 
grangian in the hierarchy. 

3. iii) In any of these hierarchies, Lagrangians depend on fields only through their 
first and second derivatives, but not on derivatives of higher order nor on the 
fields themselves. The first Lagrangian only depends on first derivatives of the 
fields. The dependence of each of the other Lagrangians on second derivatives 
is multilinear, and of order equal to the number of times an equation of motion 
has been taken to reach that level in the hierarchy. 

4. iv) All these hierarchies are finite, i.e. the iterative procedure implied by i) - 
iii) terminates after a finite number of steps. 

5. v) For each hierarchy, the last non trivial equations of motion are universal, 
namely, up to a factor, they are independent of the initial Lagrangian out of 
which the hierarchy is constructed. The associated infinite number of conserved 
charges - corresponding to the freedom in the choice of initial Lagrangian - sug- 
gests the possible integrability of these universal equations (equations of motion 
are indeed always current conservation equations for Lagrangians without an 
explicit dependence on fields). 
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Specifically, hierarchies with these properties were shown to exist in the following 
cases: 



• a single field in d dimensions, with an arbitrary initial Lagrangian (function 
of first derivatives only) [15] . 

• a single field in (d + 1) dimensions, the initial Lagrangian now being an arbi- 
trary homogeneous weight one function of its arguments 13 

• (d + 1) fields (j)"' in d dimensions, with an arbitrary reparametrisation invariant 
initial Lagrangian [15) . 

The hierarchies associated with these three cases terminate after d steps, with the 
following universal equations: 

1. 

det (f)ij = 0, (24) 

2. 
3. 

Aet{Ja4>tj) = 0. (26) 

Here, (fji and 4>ii denote the partial derivatives {d4>/dxi), {d'^(fi/{dxidxj)) of the field (j) 
with respect to the d or (d-|- 1) coordinates Xi. (The same applies of course to the fields 
and obviously the indices i and j in the equations above refer to lines and columns 
respectively of the corresponding matrices. The usual summation convention over 
repeated indices is assumed throughout). In (|26p . the quantities Ja are the Jacobians 

Ja = (~l)'^eabib2---ba(t>''l (pl^ ■ ■ ■ <!>''/■' 



Equation is a generalisation of the original two dimensional Bateman 
equation \13i (corresponding to (pSj) with (z — 1,2)) which is also known to be in- 
tegrable [13]. Finally, equation ([26]) is a generalisation to a (d — l)-dimensional 
membrane in a (d + l)-dimensional spacetime of the (universal) equation of 
motion for a parametrised particle in a flat two dimensional spacetime (corre- 
sponding to d = 1), the latter clearly being also integrable. Note that (|26p 
includes a universal equation for a string theory in three dimensions (d = 2), 
and a universal equation for a membrane theory in four dimensions (d = 3). 

Remarkably, the three classes of universal equations above are invariant un- 
der arbitrary linear GL(n) transformations in the variables Xi as well as in the 
fields (f)"', even though neither the initial nor the successive Lagrangians in the 
corresponding hierarchies would generally possess these symmetries. The equa- 
tions above thus provide examples of equations of motion admitting an infinite 
number of Lagrangians, with symmetry properties that these Lagrangians need 
not possess. 
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3.1 Other starting points 

In a second paper Hinterbickler et al [5] , the authors start from a simple massless 
Lagrangian, £2 — [tt^] in 4 dimensional space-time and generate a sequence of 
Lagrangians with second order equations of mmotion. These actually follow 
from the iterative procedure described above; and are obtained by setting each 
Ej =0, where 

f2 = -2[n] (27) 

£3 = -3([n]2 - [n^]) (28) 
£4 = -2([n]3- [n^] -3[n][n]2) (29) 

^5 = -^([n]4-6[n4]-8[n][n]3 + 3[n2]2). (30) 

The last of these, = 0, is identical, up to an irrelevant factor to equation 
([4| !.The other equations are simply given (again up to an irrelevant numerical 
factor) by 



En+i = det 



(31) 



where the matrix involved is an n x n matrix. In the above example £5 is just 
4! times the Hessian; in the previous case, £4 is just 6 times the sum of the four 
3x3 Hessians formed by the second derivatives with respect to the four possible 
subsets of three of the variables xi, X2, 2:3, X4. Likewise E^ is just 

^3 = 3( E E - a;7ra») 

p=1..4j = 1..4 

These terms are just proportional to the coefficients in the expansion of the 
characteristic polynomial of the matrix pip. These equations are solved by a 
similar technique to that above for £5 — 0. These sequences of Lagrangians and 
equations of motion are effectively the same as those of [3] . These results are 
easily understood as examples of the main results of [13] as is the fact that the 
equations of motion starting with a square root Lagrangian are of similar form 
to those starting from the Klein Gordon Lagrangian. (see [TBI). 

3.2 Multi-Galileons, 

The theory has already been extended to the case of many fields tt* [3][H] [5]- 
From the point of view presented here, just as the single-field case was a Kaluza- 
Klein reduction of the Universal Field Equation, so are the multi-Galileon fields 
a reduction of the multi-field UFE, presented in 115 , by setting the derivatives 
of the fields in the 0:5 th direction to unity. The multi-field UFE permits a first 
order formalism; The set of equations of hydrodynamic type; 

^ + E"^a^ = (32) 
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are just the equations of conservation of momentum in an inviscid, incompress- 
ible fluid. Differentiation with respect to Xj, (j = 0...n) gives a set of second 
order differential equations, and the eliminants of those together with ([5^ are 
the multi-field equations. The equations p2p admit an easy implicit solution. 
The details may be found in [TT]. Indeed, multi- field UFE , with up to n fields, 
k, say, may be obtained from the equations by setting , j — k + 1 . . .n as 
functions of u*, i . . .k, and forming the eliminants of these equations together 
with their derivatives with respect to the independent variables. In particular 
the single field UFI comes from the equation 



du ^ du ^ du _ ^ 

dxQ dxi . ^-^ dxj 



(33) 



Differentiation with respect to Xfe , fc = . . . n yields the equations 



dxo,dxk dxidxk 



j=2...n 




dxidx 



j^l...n 



Elimination of the coefiicients 1, uHu), y —u"^(u)— yields the UFI 

^-^ du f^-r 



(34) 



det 



i=l...n 

p- 

axi 

du d^u 
dxj dxidxj 



dx„ 



= 0. 



The Galileon can be obtained by a similar process of elimination, starting from 
((M)) . Making the assumption w" = 1, and also imposing the side condition 

Edu"Uu) du du^(u) 
— 1 7^ = 



m— — 1 



du dXr, 



du 



the last equation of (j34p vanishes and the eliminant is the determinant of the 
n X n matrix to zero. Thus the Galileon also has roots in fiuid dynamics. 



4 Conclusions 

The Galileon equations can be viewed as arising from a Kaluza Klein reduction 
of the universal field equation, (or more prosaically, as Chaundy has called 
it [TD], the bordered determinant) by setting the derivative of the fields with 
respect to the fifth co-ordinate to unity. This universal field equation admits 
a first order formalism in terms of an equation of hydrodynamic type, and is 
solved implicitly. The single field Galileon equation is a sum of Hessians set to 
zero and also admits implicit solutions. The universal field equation arises in 
the theory of developable surfaces, and has also connections with Dirac-Born- 
Infeld theory [17] [18] . Very recently a paper has appeared |19: showing that the 
Galileon actions and their covariant generalisations can also be obtained as a 
Kaluza Klein compactification of higher dimensional Lovelock Gravity. 
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